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In this paper, we derive two second- order of differential equation for the gluon and singlet 
distribution functions by using the Laplace transform method. We decoupled the solutions 
of the singlet and gluon distributions into the initial conditions (function and derivative of 
the function) at the virtuality Qq separately as these solutions are defined by: 

F^{x,Q^) = J-(Go,dGo) 

and 

Gix,Q^) = g{Go,dGo). 

We compared our results with the MSTW parameterization and the experimental measure¬ 
ments of {x, Q'^). 


INTRODUCTION 

The Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [1-3] evolution equations are a set of 
the integro-differential equations as can be evaluated the parton distribution functions (PDFs) 
into the evolution (where is the four momenta transfer in a deep inelastic scattering (DIS) 
process). The DIS structure functions provide our information about the deep structure of the 
hadrons and parton densities at collider and coefficient functions which contain information about 
the boson-parton interaction. These evolution equations are fundamental tools to test perturba¬ 
tive quantum chromodynamics (pQCD) in DIS experiments where the density functions can be 
evaluated by solving the DGLAP equations into the initial distributions. 

The structure function F 2 {x, Q^) reflects the momentum distributions of the partons in the nucleon. 
But the gluon distribution G{x,Q‘^) does not appear in the experimental directly, and it is deter¬ 
mined only through the quark distributions in conjunction with the DGLAP evolution equations. 
It is also important to know the gluon distribution inside a hadron at low x because gluons are 
expected to be dominate in this region. Traditionally, gluon and quark distribution functions have 
been determined using the two coupled integral-differential (DGLAP) equations. The solutions 
of the unpolarized DGLAP evolution equations decoupled have been discussed considerably over 
the last years [4-7]. There exist two main classes of approaches: those that solve the equation 
indirectly in s-space at the Laplace transform method [4-6], and those that decouple the DGLAP 
equations directly in x-space with define a hard pomeron behavior for the distribution functions [7]. 
In Refs. [4-6], authors have derived the gluon distribution and the singlet structure function into 
the initial conditions Go{x, Qq) and Fso{x, Qg). The decoupled solutions for the gluon distribution 
and the singlet structure function determined by 

Fi{x,Q^) = F{Fso,Go) and G{x,Q^) = G{Fso,Go)., (1) 

where the functions F and Q are determined by the splitting functions into the both initial con¬ 
ditions. However this method is not general, because solutions of the coupled DGLAP equations 
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are possible under certain conditions. In Ref. [7], a set of formulae to extract two second-order 
independent differential equation for the gluon and singlet distribution functions have derived. In 
this method, the singlet quark and gluon distributions have the same high-energy behavior at low- 
X according to the hard- pomeron behavior for parton distributions. Decoupling of the DGLAP 
evolution equations, for the functions F 2 {x,Q‘^) and G{x,Q‘^), determined as a function of their 
initial parameterisation at the starting scale Qo 

Fi{x,Q^) = F{dF,o,F,o) and G{x,Q^) = g{dGo,Go), (2) 

but this method is dependence to the model behavior for the parton distribution functions. 

Now we would like to find an explicit and general solution for the decoupling of the DGLAP 
evolution equations. This method is independent of the parton distribution behavior. In this 
paper, the DGLAP evolution equations coupled convert to the two second order differential 
equation which define G{x,Q‘^) and T|(x,(5^) in terms of the individual initial conditions 
{G{x,QQ),dG{x,Ql)} and {F2{x,Ql),dF2{x,Ql)} Tes])ectively. 


METHOD 


The DGLAP evolution equations can be written by 
dvr dF 2 {x,Q‘^) 


a,(Q2) 51n(g2) 


16 A - X 
4+ —In - 

6 X 




) FAx,QA\ dz 


X 


z — X 


3^ 


F|(., <3^) (l + f [ G(., (l - 2^ + 2^) f, 


( 3 ) 


and 

dvr dG{x,Q^) 1^33 —2n/ 
asiQA 91n(g2) 


-h 121n( 


1 — X, 


X 


G{x,QA + 12x I 
J X 


WG{z,QA G{x,QA\ dz 


z — x 


-|-12x f G{z,Q‘^) 

J X 


z X x^\ dz 8 

-2 + -- 


X 


z z^ 




( 4 ) 


These equations are coupled into the singlet and gluon structure functions. Let us introduce the 
variables v = ln(l/x) and w = ln(l/ 2 ;) and note that FA&~^, QA = -^1 (^) QY ^^*4 G{e~'", QA = 
G{v,QA- Using Laplace transform method in u-space [6-6,8-10], we can define the distribution 
functions in s-space as 


f{s,QA = c[Fi{v,QA;s 

and 


dw 


{w,QA', 


Fi{v,QAe-^^dv , g{s,QA = d:\G{v,QA;s\ = / Gsiv,QAe-^^dv, 

Jo 

( 5 ) 

[w,QA]\ = sg{s,QA- (6) 


= sf{s,QA , L 


dGs 

dw 



3 


We know that F 2 {v = 0, = G{v = 0, = 0 and introduce the new variable r as = 


^ for the dependence. Therefore, the coupled first order differential equations in s-space 


can be rewritten: 


and 


t) = ^f{s)f{s, t) + Qf{s)g{s, r), 


r) = ^g{s)g{s, t) + Qg{s)f{s, r), 


where the coefficient functions 0(s) and <h(s) are [4] 


(7) 




( 8 ) 


0/ (s) = 2n/ 


2 2 
+ 


+1 s+2 s+3 


(9) 


and 


, , 33-2nf fl 2 1 1 


. N - /'2 2 1 

(^) - Q 7-+ 


- + 1) - 7e , 


3V'S s + 1 s + 2 


( 10 ) 


( 11 ) 


Here ip{s) is the digamma function, 'Je = 0.5772156... is Euler’s constant and Uf is the number of 
quark flavors. Now let us find our solution for decoupling of the functions / and g from the above 
equations. The explicit form of the two second order differential equation in s-space are given by 

^(s,r) - (4>/(s)+ $g(s)) ^(s,r) - ( 0 /(s) 0 g(s) - 4 >/(s) 4 > 3 (s))r(s,r) = 0, (F = / or g). 

( 12 ) 

Solving these equations, we obtain two independent solution as 

/(s,r) = Cl exp Q(Fi(s) + i?(s))T^ + C 2 exp Q(Fi(s) - R{s))t^ , 

and 

g{s,T) = c'lexp Q(Fi(s)+ i?(s))r^ + C 2 exp Q(Fi(s) - F(s))r^ , (13) 

where Ri{s) = ^f{s) + ^g{s) , R{s) = Y^(^/(sy 77 '$^(sjj 2 Tj 74 @^^ ami the coefficients c and c' are 
the initial conditions (for the function and the derivative of function) at the Qq value. Therefore 
our solutions can be find 

/(s,r) = A:i(s,t)/(s,0) + ^ 2 ( 5 , r) |^=o. 


and 


dr 


g{s,T) = ki{s,T)g{s,0) + k2{s,T)^^^^^\r=o, 


( 14 ) 
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where 5 (s, 0 ) = ^(sjQo) /(s,0) = /(s,(5g) and the coefficient functions are defined by 


ki{s,T) = - 


and 


Ris) 


1 


(i?(s) + Ri{s)) exp ' --Ri(s)r sinh ' -R{s)t - R{s) exp ' -{R{s) + i?i(s))r 


1 


k2is,T) = 


Ris) r "V2 


2exp ( -Ri{s)t ] sinh ' -R{s)t 


(15) 

(16) 


In ti-space, we introduce two kernels Ki{v,t) and K2{v,t) with respect to the inverse Laplace 
transforms of the coefficients /ci(s,r) and A: 2 (s,r), as 

Ki{v, t) = £~^ [ki{s, t);v] , k 2 {v, r) = [k 2 {s, t);v] . (17) 

Also we can consider these two parameters at the initial scale = Qq (i.e., r = 0) where 
iLi(u,0) = S{v) and K2{v,0) = 0. Nothing that the inverse Laplace transforms of the parton 
distributions at the initial scale are given by 


Fi{v,Ql) = C-^[f{s,Ql)] , G{v,Ql) = C~^[gis,Ql)]. (18) 


Finally, we can obtain our general solutions for the decoupling of the singlet and gluon distributions 
with respect to Eqs.l4 in x-space in the following forms 


and 


G{x,Q'^)= f Ki(-,T]G{y,Ql)—+ f iLa 

^ X \?// y y X \ 2 // 


dF^{y,T)^ dy 
dr 


dG{y, t) I dy 


(19) 


( 20 ) 


or, in general form, we can rewrite the decoupled solutions as a reduced solution with the same 
coefficient functions as 


¥{x,Q^)= f Ki 

J X 


-,r)F(2/,Q2)^+ / ^2 

y J y Jx 


d¥{y,T)^ dy 


(F = Fi or G). (21) 


To get the solution this general equation, we need to the input initial conditions (Fg and clFg) 
individual for any distribution function according to the Refs. [11-12]. These results are gen¬ 
eral and gives us the exact functions of x and in a domain Xmin<x<Xmax and Q'^^n—Q'^—Qmax- 


RESULTS AND CONCLUSION 

In this method we have determined the x space distribution functions individual with respect 
to the initial and derivative of initial conditions at the stating point. We have constructed a 
general solution for the DGLAP evolution equations for the gluon and singlet distributions into 
the initial conditions respectively. To determine of the proton structure function we need to know 
the initial singlet structure functions only and also the initial gluon distribution functions for the 
gluon density which are given in Appendix A. For the validity and compatibility of the two sets of 
analytical solutions for the singlet structure function and gluon distribution function, we compared 
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the distribution functions with the HI data [14] and with the MSTW numerical solutions [13]. 
We compared the solutions of Eqs.(19) and (20) for the singlet structure function and gluon 
distribution respectively with the MSTW distributions in Figs.l and 2. Here the input for 
{F|(a;, Qq), clF|(x, Qg)} and {G(x, Qo); Qg)} are taken from the Eqs.Al-A5 corresponding 

to the lowest-Q^ point. In Eig.3 we have compared our results for the proton structure function 
with the HI collaboration [14] data at small x at different values. We consider the range 
45 < < 800 GeV^ and 0.0004 < x < 0.05 for the experimental data. In the general solutions 

of singlet structure function our computed values of F|(x,(5^) from Eq.(19) are plotted against x 
for different values of Q^. In these graphs, it is observed that the best-fit curves are obtained for 
F 2 {x,Q^) ~ as inputs for {F 2 ®(x, Qg), (9F|(x, Qg)} are taken from Eqs.Al-A3. It can 

be seen that our results are in good consistency with the HI data. 

In conclusion, we obtained two decoupling analytical evolution equation, as they are two homoge¬ 
neous second-order differential equation, for the singlet F|(x,(5^) and gluon G{x,Q‘^) distribution 
functions from the linear DGLAP evolution equations by using the Laplace transform method, 
respectively. These equations are general and require only a knowledge individual E|(x,(5g), 
Go{x,Qq) and those derivatives at the starting value Qq for the evolution, respectively. We 
observed that the general solutions are in good consistency with the experimental data and 
parameterizations. 
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Appendix A 

In Refs.[11-12] the proton structure function for 0 < x < xp has been parameterized in 
x-space as 


Ff(x,g2) 


(l-x) 


Fp 


I — Xp 


(oo ai In(g^) 02 In^(g^)) In 


xp(l — x) 
x(l — Xp) 


+ {bo + bi In(g^) -h 62 In^(g^)) In^ 


xp(l — x) 
x(l - Xp)J 


(A-1) 


where xp = 0.09 and F^ixpjQ'^) = Fp = 0.41. Table I shows fitted quantities and their errors. In 
region x < 0.1 we can approximately write F|(x,g^) ~ IS/hiT^(x,Q^). So in Qq = 1 we have 


dFS , ^ 2 m nT 1 f xp{l — x)\ 2 /^ 2 :p(l —x)\ 


or 


9^2, 2m nT , f Xp{l — x)\ 2 f Xpi'i- 

-^(x. Q c ^^(1 - a:) [ai In (j + 6. In 


— x) 


- xp) 

Also, the LO gluon distribution function G{x,Q‘^) = xg{x,Q‘^) has been determined by 


(A-2) 


(A-3) 


G(x, Q^) = cq + ci In(g^) C 2 In^(g^) + [do + di In(g^) d-2 In^(g^)] In 


X 


[eo + ei ln(g2) -b 62 In^(g^)] In^ . (A-4) 

The parameters have been shown in Table I. The derivative of LO gluon distribution function in 
Qq = 1 becomes 

g(==.<?blo..Q2 = ^ (ci + rfl In (1) + e, In^ (1)) , (A-5) 



TABLE I: The fitted parameters 


parameters 

value 

parameters 

value 

Co 

-0.459 

ao 

-5.381 X 10-2 ±2.17 X 10-3 

Cl 

-0.143 

ai 

2.034 X 10-2 ± 1.19 X 10-3 

C2 

-0.0155 

02 

4.999 X 10-3 ±2.23 x 10-* 

do 

0.231 

bo 

9.955 X 10-3 ±3.09 x 10-* 

di 

0.00971 

bi 

3.810 X 10-3 ±1.73 X 10-* 

d2 

-0.0147 

b2 

9.923 X 10-4 ±2.85 x 10-^ 

Co 

0.0836 



Cl 

0.06328 



62 

0.0112 
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FIG. 1: The singlet structure function compared with the MSTW parameterization [13] (solid line) at = 
20,100 and Ml GeV‘^( bottom to top, respectively). 
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FIG. 3: Proton structure function Q^) as function of x at different values compared with the HI 
data [14]. 



